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Abstract-A compact template for the division of two giving 
polynomials to find quotient and reminder is derived. The 
process is very simple and direct, comparing to the familiar 
classical long polynomial division and synthetic polynomial 
division. 
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I. Introduction 

There are several approaches for finding the quotient and 
remainder from dividing two given univariate polynomials. 
Long polynomial division is very popular but tedious in 
computation, and widely used even by high school students. 
Synthetic polynomial division is fairly easy to use but only 
appropriate for the linear divisor [1]. Convolution polynomial 
division [2] is direct in operation, and used in MATLAB 
built-in routine. 

This work presents a compact template for polynomial 
division. The process is very simple and straightforward and 
does not need to write down any intermediate steps, as in the 
familiar classical long polynomial division and synthetic 
polynomial division. It is extremely suitable for hand 
computation with a plain calculator. 

The concept of compact template may be easily applied to 
the product of two given univariate polynomials. 

Two typical numerical examples are provided to show the 
merit of the approach presented 

II. POLYNOMIAL DIVISION 

The division of two given polynomials is expressed as 

a(x) a(x) 

where the dividend b(x) of degree n and divisor a(x) of degree 
m are 

b(x) = b x" + b^"- 1 + + b n _ t x + b n 

a(x) = a x m + a 1 x m - 1 + + a m _ x x + a m 

Then the quotient q(x) of degree n-m and remainder r(x) 
of degree m-1 are obtained, 

q(x) = q x- m + qy- 1 + + q n _ m _,x + q n _ m 

r(x) = r Q x m -' + r lX m - 2 + + r m _ 2 x + r m _, 

By carefully examine the familiar classical long 
polynomial division we may find the following relations 
among the coefficients of these polynomials: 

K = Qo a k + Qi a k-i + + Q k -i a i + Qk a o> 

k = 0, !,•••, n-m 



b k -r k _ (n _ m+1) =q Q a k +q 1 a k _ 1 + + q n _ m a k _ (n _ m) , 

k = n-m + 1, • • • , n 

The computation process can then be made by the 
following "compact polynomial division template" 
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and the desired coefficients are therefore calculated from the 
given coefficients in this template, 

q k =([l]-[bj -[q .-'-.q k -i]-[ a k .---. a i]) / a o. 

k = 0, 1, • • • , n - m 
r k =[l]-[b nMk ]-[q ,---,q n - m ]-[a n _ m+1+k r-,a 1+k ], 

k =0,1, — , m-1 

Here [ ] • [ ] denotes the inner product of two equal 
size arrays; and it is clearly understood that 
b l = 0, I > n, and a, = 0, I > m. 

The quick check can readily be achieved by the following 
simple operations: 

(a + a 1 + --- + a m )(q + q 1 + --- + q n _J 
= (b + b 1 + --- + fa n )-(r + r 1 + --- + r m _ 1 ) 

If all the given coefficients are integers and d Q ^ 1 , the 

scalar multiplier 1 in the template may be replaced by a^ (n- 
m+1), the desired coefficients are then obtained in pure 
arithmetic operations such that the round errors may be 
eliminated. 

The compact polynomial division template for 
n = 8, m = 5 is shown for illustration: 



"0 



h h h 
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j % a l a 2 Q 3 Q 4 °5 

"1 ^ q\ q~ 2 % | r Q \ \ \ r~ 

The desired coefficients are computed in order: 

q = ([l]'[b D la 

q 1 = ([UTbJ - [q ][ Q i]) /a o 

q 2 = ( [ 1 ] ■ [b 2 ] - [ q q 1 ] • [ a 2 a 1 ] ) / a 

q 3 = ( [ 1 ] • [f> 3 ] - [q q t q 2 ]-[a 3 a 2 a 1 ]) / a 
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[l]-[i> 4 ] - [q q t q 2 q 3 ] • [ a 4 a 3 a 2 a i 1 
r t = [ 1 ] • [b 5 ] - [ q q 1 q 2 q 3 ] ■ [ a 5 a 4 a 3 a 2 ] 
r 2 = [ 1 ] ■ [b 6 ] - [ q a q 1 q 2 q 3 ] ■ [ a 5 a 4 a 3 ] 
r 3 = [ 1 ] • [b 7 ] - [q q 1 q 2 q 3 ][ a 5 a 4 ] 
r 4 = [ 1 ] ■ [b 8 ] - [ q q x q 2 q 3 ] • [ o 5 ] 

III. NUMERICAL EXAMPLES 

Example!. The division of the two given polynomials: 
b(x) = 4x 8 + 5x 7 - x 6 + 7x 5 - 6x 4 + x 3 + 2x 2 - 3x + 7 



a(x) = 3x 5 + x 4 - 7x 3 + 5x 2 - 4x +2 



will gives the quotient and remainder: 



q(x) 
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r(x) = x H x x 
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The desired results may be determined either one of the 
following approaches: 

From the comparison of the approaches in this example, 
the convolution polynomial division is obviously much simple 
and effective. The desired quotient and reminder are readily 
determined without computing any intermediate values as in 
the familiar classical longhand polynomial division and 
synthetic polynomial division. 

l.Long polynomial division: 

2. Synthetic polynomial division 

3. Compact polynomial division template: 

4. Compact template in pure integer operation: 

(l)Long polynomial division: 
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2)Synthetic polynomial division 
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(3)Compact polynomial division template: 
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(4)Compact template in pure integer operation: 

44 +5 -1 +7 -6 +1 +2 -3 +7 
4-) +3 +1 -7 +5 -4 +2 
81 +108 +99 +192 +176 +619 +533 -143 +77 +215 

Example2. A given polynomial is to be expressed as 
powers of linear factor (2x -3): 



p(x) = 2x 7 -13x 6 + 51x 5 -119x 4 +191x 3 -181x 2 + 104x +6 



+ J-(2x-3) 7 + i(2x-3) 6 +^(2x-3) 5 +-^(2x-3) 4 
64 8 64 16 



^ 811 ^ ^3 125,- ,,,3259... ,,1125 

+ — (2x-3) +— -(2x-3) + — — (2x-3) + — — 

64 4 64 16 

by working the polynomial divisions consecutively: 
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IV. POLYNOMIAL MULTIPLICATION 

The concept of compact template may be easily applied to 
the polynomial multiplication. Let a(x) of degree m and b(x) 
of degree n be the given polynomials defined as before. The 
resulted polynomial p(x) of degree m+n is the multiplication 
of these two given polynomials: 

p(x) = a(x) ■ b(x) 

or (P ( / M +P/ W " 1 +-'-+P„ M -i x +P™ ) 

= (V+a l x , ^+---+a^x +a m )-(b/ +b 1 x" 1 +---+b_ 1 x +b ) 

The desired coefficients are then determined from the 
given coefficients, 

Pk = a oK + a A-i + ••• + a k -A + a A> 

k = 0,1, ••• , m + n 

The computation process can be easily made by the 
following "compact polynomial multiplication template" 



and the desired coefficients are thus calculated directly from 
the given coefficients in this template, 



p = [a , ■ ■ ■, a 1 • [b , • • •, b 1, k = 0, 1 ■ • •, m + n 

1 k k k a 

Here [ ] ■ [ ] denotes the inner product of two equal 
size arrays; and it is clearly understood that 
b ( = 0, I > n, and a j = 0, I > m. 

V. CONCLUSION 

The useful template is derived for division of polynomials. 
By comparison from the examples, this approach is simple 
and effective. The desired quotient and reminder are readily 
determined without writing down any intermediate values as 
in the familiar classical longhand polynomial division and 
synthetic polynomial division. 

One of the important applications is to find the roots with 
multiplicities of any given polynomial after the GCD of the 
polynomial and its derivative is computed [3-4]. 
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